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We implement our past investigations in the quark-antiquark interaction through a 
non-perturbative running coupling defined in terms of a gluon mass function, similar 
to that used in some Schwinger-Dyson approaches. This coupling leads to a quark- 
antiquark potential, which satisfies not only asymptotic freedom but also describes 
linear confinement correctly. From this potential, we calculate the bottomonium and 
charmonium spectra below the first open flavor meson-meson thresholds and show 
that for a small range of values of the free parameter determining the gluon mass 
function an excellent agreement with data is attained. 
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I. INTRODUCTION 


The use of perturbation theory in quantum chromodynamics (QCD) [I| provides a useful 
tool to study high energy processes as conhrmed by the accurate description of deep inelastic 
lepton-nucleon and nucleon-nucleon collision processes j2^- 

Asymptotic freedom 0,0 is the fundamental property of QCD that justihes the use of 
perturbation methods to describe the interaction among the constituents at high momentum 
transfers (short distances). 

In order to describe low energy (large distance) physics, for example the hadron spec¬ 
trum, perturbation theory cannot be used. The reason being that the interaction among 
the constituents is dominated by another fundamental property of QCD, the confinement 
of quarks and gluons 0], which is highly non-perturbative. Efforts have been made and 
techniques developed to study this low momentum transfer regime, e.g. lattice QCD js J9|, 
phenomenological quark potential models [l^ 12|. studies on the nontrivial vacuum 13], 


effective theories dispersive extensions of the beta function [15|, [1^, and the solution 
of the Schwinger-Dyson equations 00. Another development which has led to very suc- 
cesful results in the study of the spectrum and the behavior of the strong coupling in the 
infrared (IR) is AdS/CFT 0,0. 

The resolution of Schwinger-Dyson equations (SDE) leads to the freezing of the QCD 
running coupling (effective charge) in the infrared which can be conveniently understood as a 
dynamical generation of a gluon mass function, giving rise to a momentum dependence which 
is free from IR divergences 0, 18, 21-23] . It was shown that the interquark static potential 
for heavy mesons described by a massive one gluon exchange interaction obtained from the 
propagator of the truncated Schwinger-Dyson equations does not reproduce the expected 
Cornell potential HSl. This was attributed to the lack of a mechanisms for conhnement 


by explicit comparison with lattice QCD calculations [25 


To solve this problem an ad hoc effective mechanism was proposed based on a singular 
nonperturbative coupling 0], which generates a Gribov singularity in the potential while 
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keeping the propagator infrared finite. We show here that the proposal leads to a potential 
with linear confinement and the correct asymptotic behavior of QCD. The resulting potential 
is Cornell like for intermediate distances (0.1 — 4 fm) and for a small range of values of the 
free parameter determining the gluon mass function it provides a good spectral description 
of heavy quarkonia. 

Other authors have also used non-perturbative effective couplings which reproduce 
asymptotic freedom in the UV and have 1/q^ behavior in the IR to study the spectrum and 
other physical observables 27-3^1. These couplings have been justified from QCD under cer- 

However, it must be noted that many succesful phenomenological 


tain requirements |33H35 


studies of the IR coupling point to a finite value 

The contents of this paper are organized as follows. In Section Ull we establish the pre¬ 
scription to get the static quark-antiquark potential from the running coupling. A possibility 
to generate a linear confining potential is through a strong Gribov singularity for small mo¬ 
menta. We propose a functional form resulting from a modification of the Schwinger-Dyson 
solution for the coupling, which has the desired singularity structure at small momenta and 
satisfies asymptotic freedom at large momenta. The presence of the singularity requires, in 
order to get the potential, a regularization scheme which is described in detail in Section 
uni Then, in Section HVl the results obtained are presented. The potential is compared to 
the Cornell one. The calculated heavy quarkonia masses are compared with data and with 
the ones resulting from “equivalent” Cornell potentials. Finally, in Section |Vl the main 
conclusions of this work are summarized. 


II. QUARK-ANTIQUARK POTENTIAL FROM THE NON-PERTURBATIVE 

RUNNING COUPLING 


The history of the theoretical description of the quark-antiquark potential and its relation 
to QCD has been motivated by the tremendous phenomenological success of the Cornell 
potential [13, 11 

r) = -- + ar , (1) 


U 


c. 


r 

where y, the Coulomb strength, and a, the string tension, are constants to be fitted from 
data. It is amazing that this simple potential reproduces quite accurately the experimen¬ 
tal heavy quark meson spectra below the open flavor meson-meson threshold energies (see 
Section ED. 

Two limits characterize this potential. The short distance limit which satisfies asymptotic 
freedom, and the long distance limit which describes the IR behavior, i.e. confinement of 
heavy static sources. 

In general, a quark-antiquark potential in configuration space, V{r), reads (note that the 
angular variables have been integrated) 


V{r) = — / V{q^)q 


2 2sin(gr) 


TT 


dq , 


'0 Qr 

where V{q^) is the potential in momentum space. 

By assuming one gluon exchange dominance, V{q^) can be derived from QCD as 


( 2 ) 


V{q^) = -«.(g^)A(g2) 


( 3 ) 
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with Q!s(g^) and A(g^) standing for the running coupling and the gluon propagator respec¬ 
tively. 

It is now well established from lattice QCD, and confirmed by Schwinger-Dyson calcula¬ 
tions, that the finite character of the propagator (at 0) reads 00 


A(g^) 


1 

-|- w?g{q^) 


(4) 


where mg{q^) is an effective gluon mass (see below). Other points of view to parametrize 
the gluon propagator have been developed, like flux tubes 0] and effective potentials j0 . 
in which the connection to the concept of massive gluon exchange is not currently known. 

Regarding the coupling, a truncated solution of a gauge invariant subset of the Schwinger- 
Dyson equations for QCD gives rise to a functional form showing a freezing in the IR which 
can be cast as 0] : 


a 


(SD) 


(a = 


47r 


/doln 


q2+m2(g2)\ 


(5) 


A2 


where I3q = 11 — 2/3 nj is the first /9-function coefficient for QCD, rif being the number of 
active quarks, A is the scale parameter in QCD and m?g{q^) is a gluon mass given by 


ml{q^) = 


rrin 


1 -h 


( 6 ) 


with mo = mg{q^ = 0), Ad and p > 0 constants. Other parametrizations for both coupling 


and gluon mass can be found in ref. [36[. The coupling (|5]) is to first order the one proposed in 
"HI obtained from the gluon propagator using the pinch technique. Note that the dynamical 


23| and therefore is determined by 


gluon mass varies with the number of flavors {rif) |! 
the initial mo and the QCD scale A. 

The combination (I3j), (0]), (|5]) and ([6]) is renormalization group invariant since the ghost 
contribution has been incorporated into the definition of the coupling [24. Therefore one 
can construct a potential. The result is a non linearly rising potential [24 contrary to the 
expectation from the quenched approximation followed to derive the propagator and the 
coupling. Actually a linearly rising Cornell-like behavior has been obtained in quenched 
lattice QCD j§, 25|. One possible reason for this anomaly might be that confinement is 
more than a one gluon exchange effect. An alternative possibility, suggested in reference 


26 


is that the Schwinger-Dyson coupling contains the physics at intermediate and large g^ 
but it lacks some vertex corrections at low g^. We shall examine this alternative next. 

For this purpose we realize that in order to lead to linear confinement the mass singularity 
of the propagator must be eliminated and instead a Gribov singularity should appear. An 
easy way to do this from (|5]) is through the modified coupling 


as{qy 




(g^ + m2(g2)) 


(7) 


so that the factor (g^ -|- ni^^qy'j cancels the singularity of the propagator whereas the factor 
gives rise altogether with (ai®°^(g^)) to a A Gribov singularity. In this way the 
low g^ behavior of a^^^\qy has been corrected as required while the intermediate and long 














4 


distance dependencies are preserved since is a qnickly decreasing function with 


When mo = A our potential becomes in the IR limit, << A^, that of Richardson 27 


signaling that it is dehned in the MS scheme. Moreover, Q!s(g^) reproduces in the asymptotic 
limit g^ >> A^ the well known one-loop perturbative QCD result 


arV) = 47r/(Alti(,7A=)) , 

Additional terms of the order C>(l/(ln(g^/A^)g^)), where N = A + 2p > A, have to be 
added to this result but hey are negiglible and arise from the IR behavior js^. 

In order to check whether this modihed coupling proposal is physically meaningful or not 
we shall study, from the potential deriving from it, heavy quarkonia, using as a criterion 
of meaningfulness the accurate description of the spectra. More precisely, we shall choose 
for the evolution mass parameters Ai and p their estimated values in the Schwinger Dyson 
context 

M. = 436 MeV and p = 0.15, 

while leaving A (= mo) as the only free parameter of the potential. Then, an accurate 
spectral description for a reasonable value of A might be an indication that the proposed 
coupling parametrizes efficiently and conveniently the phenomenology. 

It must be stressed at this point that we have incorporated the Gribov singularity guided 
by lattice QCD and simplicity and in so doing we have departed from the SDE approach. 
Nevertheless the effective way of introducing the coupling mantains the behavior of the SDE 
result for large and intermediate energies and therefore a good ht of the spectrum is strong 
support of the SDE calculation. 


III. REGULARIZATION PROCEDURE 


From ([2]), dl]) and ([7]) the potential reads 


U(r) 


32 

Wo 


sin(gr) 


In 


„2_i_ 

^ l + {q'^ / 

A2 


qr 


dq 


( 8 ) 


The integral in ([8]) is divergent since the integrand behaves as ^ at the origin g 0. To 
extract its physical content we have to regularize it. To understand how the regularization 
procedure works let us consider the simpler case of the potential (with the same singular 
behavior at the origin) 

4z/A2 

where u is an dimensionless constant. Its Fourier transform, from ([2]), is given by 


U(g2 


V^{r) 


8z/A^ 1 sin(gr) 
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By introducing a regulator A we can rewrite it as 


o 1 

V^{r) =-;-lim 


Svr r A^o Jg _j_ 
8 z/A 2 1 


Stt r A^-o 


lim ( Im 


sm.{qr)dq 
Q 


1-00 (q^ + A2)2 V 


so that the integral can be solved analytically in the complex plane giving 


V^{r) = 


2 z/A 2 


lim ( — 

A-i-O 


A 


2 z/A 2 


The divergent term limA^.o (~y) does not depend on r. As the potential is dehned up 
to an arbitrary constant we may simply remove such spurious divergence from the potential. 
To do this in a more systematic way, let us realize that the divergent term comes out from 
the behavior of the integrand at 0. Then, by making an expansion of the integrand 

around g = 0 , 

1 sin(gr) 1 (gr)^ 

g 2 q2 Qq2 ’ 

we can easily identify the hrst term ^ as causing the divergence of the integral at g^ 0 . 
Indeed, if we integrate this term with the same regularization procedure we get 




jdq = 


2 z/A 2 


3vr Jo (g2 + A2)' " 3A 
Therefore we can redehne the potential by subtracting this non physical divergence as 


V(r) ^ -hh! lim 


sin(gr^ 


Stt a^o \Jq (g 2 x2y 


qr 


dq 


lo (g^ + A 2 ) 


rdg 


Back to ([H]) it is more convenient to regularize it through a cutoff 7 in the form 


V{r) = 


4 lim 

3/^0 7-^0, 


sin(gr^ 


In 


A 2 


^ 1 + (92/A42)1+P 

A 2 


qr 


dq 


To subtract the spurious divergence we expand the integrand around g = 0 : 

1 sin(gr) 


In 


„2_i_ a2 

^ ~^1 + ( 92 / 7 V | 2 ) 1 + p ' 


A 2 


A" , 


A 2 


Af 2 + 2 P 


qr 


g'^ + 


A^ 




g'^ + 


A 6 




g®^ +... 


and keep only the terms giving rise after integration to a singular behavior at g^ 0 (they 
correspond, for p = 0.15, to the explicitly written ones). By integrating these terms with 
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FIG. 1: Bottomonium (solid) and charmonium (dashed) potentials for A = 320 MeV. 



Kfm) 


the chosen cutoff regularization we get 


-q^^ + 




-q^P + 


A 6 


Is (7) _ r 1 A , 

A 2 q^\ ' M^+^p'" ' M^+^p 

1 A2 1 1 A^ 


q^P ] dq 
1 1 


^ _A^2+2p _ 2p) ^l-2p ;\/f4+4p ^l-4p 

A® 1 1 


_A/f6+6p (1 - 6p) 7^-6^ 
so that the physical potential reads 

/ 

32 

V(r) = -— lini 

3/9o 'y~^^ 


\ 


sm{qr) 


In 


\ qr 

A 2 


dq - Is ( 7 ) 


/ 


(9) 


which is evaluated numerically. 


IV. RESULTS 

In order to £x A, the only free parameter of the potential, we require (jO]) to provide 
a reasonable description of the heavy quarkonia (bottomonium and charmonium) spectra. 
As will be justihed later on, for bottomonium, one has to use n/ = 4 and for charmonium 
Hf = 3. To calculate the spectrum we solve the Schrodinger equation. For any value of A, 
we choose the quark masses, mb and me, to get the best spectral £t. In this regard, as ([9]) 
represents a quenched potential we restrict the comparison with data to energies below the 
corresponding open flavor meson-meson thresholds. 

It turns out that only for a quite restricted range of values of A (Ac ~ A;, ~ 320 MeV) a 
good spectral description for bottomonium and charmonium is obtained. The corresponding 
potentials for A = 320 MeV are shown in Fig. [1] 
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jPC ^%(r)A=320 MeV ^^PDG 

MeV MeV MeV 


f— Is 

9489 

9460.30 ± 0.26 

9479 

2s 

10023 

10023.26 ±0.31 

10013 

Id 

10147 

10163.7 ± 1.4 

10155 

3s 

10354 

10355.2 ±0.5 

10339 

2d 

10435 


10427 

4s 

10621 

10579.4 ± 1.2 

10596 

3d 

10681 


10666 

5s 

10854 

10876 ±11 

10825 

4d 

10903 


10883 


(0,1,2)++ 

Ip 

9903 

9899.87 ±0.28 ±0.31 

9920 

(0,1,2)++ 

2p 

10254 

10260.24 ± 0.24 ± 0.50 

10252 

(0,1,2)++ 

3p 

10531 

10534 ± 9 

10519 


TABLE I: Calculated bottomonium masses from V (r) A =.‘H?n MeV and rrib = 4450 MeV. Masses 
for experimental resonances, Mpdq, have been taken from 39|. For Ip and 2p states the exper¬ 
imental centroids are quoted. For 3p states the only known experimental mass is listed. Masses 
from {r)^^ , the “equivalent” Cornell potential (llOp and the same quark mass are also shown 
for comparison. 


As can be seen, the potential ([9]) shows a soft flavor dependence in the slope for interme¬ 
diate and large distances. However, if we use n/ = 3 also for bottomonium the calculated 
spectral masses will only change slightly. 

In Tables HI and UTI we list the calculated masses for bottomonium and charmonium for A = 
320 MeV as compared to data. To denote the states we use the spectroscopic notation nl, 
in terms of the radial, n, and orbital angular momentum, I, quantum numbers of the quark- 
antiquark system. As we are dealing with a spin-independent potential we compare as usual 
the calculated s—wave state masses with spin-triplet data, the p—wave state masses with 
the centroids obtained from data and the d—wave states with the few existing experimental 
candidates. 

A very good spectral description is attained. Note that the only signihcant difference 
(> 60 MeV) between the calculated masses and data is for the 3s charmonium state and it 
may be explained through conhguration mixing with the 2d one. 

For further comparison the spectra from “equivalent” Cornell potentials have also been 
quoted. This equivalence is based on the observation that for intermediate distances (0.1 — 4 





nl 

Mv(r)A=320 MeV 

MpDG 

MyCp)^ 

MeV 

MeV 

MeV 

1— 

Is 

3090 

3096.916 ±0.011 

3105 


2s 

3671 

3686.108l[];[]}^ 

3671 


Id 

3772 

3778.1 ± 1.2 

3772 


3s 

4102 

4039 ± 1 

4096 


2d 

4170 

4191 ±5 

4166 


(0,1,2)++ Ip 

3484 

3525.30 ±0.11 

3493 

2++ 2p 

3940 

3927.2 ±2.6 

3940 


TABLE II: Calculated charmonium masses from V (r) A= 33 n MeV and rric = 1030 MeV. Masses 


for experimental resonances, Mp£)g, have been taken from [39|. For Ip states the experimental 
centroid is quoted. For 2p states we quote the state which lies below the threshold. 
Masses from V^{r)cc , the “equivalent” Cornell potential (1121) . and the same quark mass are also 
shown for comparison. 


fill) the potentials ([9]) can be very well approximated by the Cornell types 

= + ( 10 ) 

with 

= 800 MeV.fm~\ (11) 

Xbi = 100 MeV.fm, 
a^l = 693 MeV, 

and 

{V^{r))^_ = accr-^ + acc ( 12 ) 

with 

<Jcc = 735 MeV.fm“\ (13) 

Xcc = 100 MeV.fm, 

Occ = 658 MeV, 


where the htted Coulomb strengths are in agreement with the values derived from QCD 
from the hyperhne splitting of Ip states in bottomonium 40| and from the hne structure 
splitting of Ip states in charmonium j^ . 

As a matter of fact, in the mentioned region, they can not be distinguished from those 
in Fig. [U However, below and above this region they become different as shown in Figs. [2] 
and [3] for bottomonium. 
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FIG. 2: Bottomonium potential for A = 320 MeV (solid line) versus its “equivalent” Cornell 
potential (dotted line) for short distances. 



FIG. 3: Bottomonium potential for A = 320 MeV (solid line) versus its “equivalent” Cornell 
potential (dotted line) for long distances. 


These differences have not much effect on the calculated masses, as can be seen in Tables [T] 
and m However, in order to get a more accurate £t to the spectra from a Cornell potential 
the string tensions have to be slightly increased with respect to those of flTOll and (El)- 
It is interesting to check a posteriori our initial assumption about the number of active 
quarks rif {3 for charmonium and 4 for bottomonium). The momentum determines the 
active number of flavors in the coupling, thus if 

rrifC <q ={Pq-Pq) =4|pg| 

where is the three-momentum of the quark (charm or bottom) in the center-of-mass 
system, then rij is to be used in the coupling. This can be conveniently rewritten as 

2 2 , 22 
rrijC < m^v 
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Kfm) 

FIG. 4: Charmonium potential (solid line) for A = 320 MeV (n/ = 3) as compared to Richardson’s 
potential (dashed line) for A = 398 MeV (n/ = 3). 


where v is the relative quark-antiquark velocity. By using the values calculated from the 
Is wave functions, = 0.3c^ and = 0.09c^, we can see that in charmonium 

q > mgC, s for strangeness, and therefore rif = 3, while in bottomonium, q > rricC and 
therefore rij- = 4. 

It is illustrative to compare the potential (|9]) htting the charmonium spectrum (with 


A = 320 MeV, see Table [IT)) with the potential of Richardson 27 
A(g2) = 4 and 


obtained from 


a 


(R) 




dvr 


with 


(14) 


/3oln(l -FgVA2) 

This comparison is shown in Fig. IH Although generated from different approaches, 
both potentials give a similar quality hts to the spectrum. It must be recalled that both 
potentials are dehned in the MS scheme. The value of A required by equation (Ej) (A = 320 
MeV) is in better agreement with QCD than Richardson’s (A = 398 MeV) [1^, since in the 
MS scheme, A (4 Loops, rif = 3) = 358 ± 22 MeV and A (4 Loops, rij = 4) = 303 ± 21 
MeV Ei-El. 


V. CONCLUSIONS 


This work has been motivated by previous studies aiming at describing the phenomeno- 

Our aim here has been to describe the heavy 


logical successful potentials from QCD |24126 


quark spectra from a potential derivable from non-perturbative QCD studies. Our research 
has led to a simple description of linear conhnement in quenched QCD in terms of a gluon 
mass function describing the non-perturbative coupling. 

There are several issues which have arisen in our investigation that merit attention and 
which we next recall. 

Conhnement can be described by a one-gluon-exchange picture if some of the long¬ 
distance physics is folded in an effectively generated gluon mass. In that context asymptotic 
freedom of the coupling and phenomenology require that the gluon mass function goes 
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rapidly to zero, faster than 1/q^. This result is in complete agreement with lattice QCD and 
the SDE solutions. 

In that context, a Gribov type singularity is the least requirement for linear conhnement. 
This IR singular behavior is an ad hoc assumption, which many phenomenological studies 
of the IR coupling do not support [s^, and implies that the SDE formulation of the gluon 
mass function must be restricted to mo = A. If mo > A, linear conhnement will soften 
to a Yukawa type behavior. In the present non-relativistic dynamical scheme this type of 
potentials are non conhning, however this is not so in other formulations j^ . 

A main result of our analysis is that we are led to Cornell type potentials in the region 
relevant to describe the spectra (0.1 — 4 fm). In this way we support the phenomenological 
success associated with Cornell type potential as a consequence of QCD. 

Our £t of the spectra is excellent. Note that we have used only one free parameter, 
which comes out at a very reasonable value (A ~ 320 MeV). Moreover, we justify the slight 
difference between the bottomonium and charmonium potentials in terms of the number of 
havors entering the coupling. 

Therefore we have shown that by incorporating the infrared Gribov singularity in a man¬ 
ner that respects the behavior of the massive SDE coupling with hxed parameters at large 

we obtain an excellent potential capable of reproducing the heavy hadron spectrum with 
only one parameter. A, which moreover comes at a wishful value ~ 300 GeV close to the 
value determined from other phenomenology in the MS scheme. Moreover, the havor de¬ 
pendence of the SDE coupling controls the havor dependence of the potential. While our 
potential is very close to the Cornell potential in the physical range it dihers from it outside 
the range and has a perfect asymptotic QCD behavior. Moreover this coupling not only 
can be used to explain the spectrum but can be used in many other processes where we are 
dealing with large to moderate energies. 

We conclude by stating that we have found an explanation for the phenomenological 
successful potentials in terms of a non-perturbative ehective mechanism describing the strong 
coupling. This mechanism is defined by means of a gluon mass function with properties 
closely related to lattice QCD and SDE studies to which we have incorporated in a gentle 
way the Gribov singularity. We have supported the successful Cornell type potentials as 
arising from specihc mechanisms in non-perturbative QCD. 
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